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I. INTRODUCTION 

The preparation of ultracold atomic Bose and Fermi gases 
in various trapping environments allows to study in a pre- 
cise way important aspects of quantum many-body dynam- 
ics |[ll|2l[3, 4, 5, 6, 7, 8]. For this reason the field has attracted 
in recent years researchers from a variety of physical disci- 
plines, ranging from condensed-matter to high-energy particle 
physics and even cosmology. In past experiments with Bose- 
Einstein condensates of dilute gases it has been found that 
these are, in many cases, approximately described by a com- 
plex scalar field which solves the Gross-Pitaevskii equation 
(GPE) f9', T^. Despite the fact that the first-order coherence 
reflected by this equation has its origin in the quantum nature 
of the Bose condensation phenomenon, the GPE arises as the 
classical field-theory approximation of the underlying quan- 
tum many-body problem. It thus neglects all quantum sta- 
tistical fluctuations contributing to the dynamics of the scalar 
field. However, it is the role of these quantum statistical fluc- 
tuations which is of central importance for our quantitative 
understanding of a wealth of phenomena described by quan- 
tum many-body dynamics. Accordingly, experiments which 
are sensitive to fluctuations are crucial to test our theoretical 
understanding of complex many-body problems. 

If fluctuations are relevant then the quantitative interpreta- 
tion of the data typically requires nonperturbative theoretical 
descriptions, which often have to take into account nonequi- 
librium dynamics to match realistic experimental situations. 
Two cases should be distinguished in this context: If the real- 
time dynamics of a Bose gas is dominated by classical statis- 
tical fluctuations then it can be well approximated by a large 
number of numerical integrations of the classical field equa- 
tion (GPE) and Monte Carlo sampling techniques lfTll[T2l[l3ll . 
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This takes into account nonperturbative dynamics, however, 
neglects all quantum corrections. A corresponding classical 
statistical description does not exist for fermions. The other 
case concerns dynamics where quantum fluctuations are rele- 
vant. It is an important challenge to quantitatively determine 
the role of quantum fluctuations, thus predicting signatures 
for the detection of genuine quantum effects in an experiment. 
Here, we address this question for a Bose gas. 

In fact, many experiments concerning ultracold Bose gases 
fall short of being sensitive to quantum statistical fluctuations, 
and can be accurately described by the GPE. The importance 
of classical statistical fluctuations can rise if the gas is suffi- 
ciently dense. A combination of low densities and strong self- 
interaction can lead to enhanced quantum fluctuations as com- 
pared to classical statistical fluctuations. Zero-energy scatter- 
ing resonances, particularly the so called magnetic Feshbach 
resonances Q III [M O H so far have played a leading 
role in the creation of strong interactions in degenerate atomic 
quantum gases. Near a Feshbach resonance, the scattering of, 
e.g., a pair of Bose-condensed atoms, whose relative energy 
is very close to zero, can be described by a strongly enhanced 
s-wave scattering length a. Present-day experimental tech- 
niques allow for resonance-enhanced scattering lengths larger 
than the mean interatomic distance {N/V)^^^"^ in the gas. As 
a consequence, the diluteness parameter a'^N/V is larger than 
one. The Bose-Einstein condensate is no longer in the colli- 
sionless regime, it represents a strongly interacting system. 
Feshbach resonances have gained a strong practical impor- 
tance for fermionic gases, where losses are suppressed in the 
unitary limit [19.1 and where they allow to study the transition 
from a phase of Bose-condensed molecules to a BCS type su- 
perfluidSlill^lIll]. 

One- and two-dimensional traps ll22l l23ll as well as opti- 
cal lattices ll24ll25ll allow to reaUze strongly correlated many- 
body states of atoms. In an optical lattice, strong effective 
interactions can be induced by suppressing the hopping be- 
tween adjacent lattice sites and thus increasin g th e weight of 
the interaction relative to the kinetic energy l24l l26ll . This 
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leads, in the limit of near-zero hopping or strong interactions, 
to a Mott-insulating state |27]. In a one-dimensional trap, the 
gas enters the so-called Tonks-Girardeau regime, if the dimen- 
sionless interaction parameter 7 = gi-[:imV/{h'^N) is much 
larger than one |28, 29, 30], Here, t^iD is the coupling param- 
eter of the one-dimensional gas, e.g., giu — 2Y?al (mlY) for 
a cylindrical trap with transverse harmonic oscillator length 
l^. In the Tonks-Girardeau limit 7 — > 00 the atoms can no 
longer pass each other and behave in many respects like a one- 
dimensional ideal Fermi gas |23J. 

The theory of the full nonperturbative real-time quantum 
dynamics is in general a demanding problem. Akeady the 
description of weakly correlated many-body dynamics suf- 
fers from the problem that it requires summations of infinite 
series of perturbative processes. These summations can be 
efficiently taken into account using functional-integral tech- 
niques for the quantum field theory, which are based on the 
two-particle irreducible (2PI) effective action ifsil l32i [33ll . 
Much progress has recently been achieved using nonperturba- 
tive expansions of the 2PI effective action to next-to-leading 
order in the number of field components |34, 35]. This has 
been employ ed in the context of ultracold quantum gases in 
Ref. ll36i Wk Issll . It has previously been successfully used 
to study far-from-equilibrium dynamics and thermalization in 
relativistic bosonic |34, 39, 40, 41] and fermionic ]42, 43] the- 
ories. For an introductory review see Ref. In Ref. ^ 
the approach has been used to compare quantum and classi- 
cal statistical nonequilibrium dynamics for a relativistic scalar 
field theory in the absence of a field expectation value. Since 
the dynamic equations for the quantum and classical corre- 
lation functions differ only by few characteristic terms, this 
can be used to derive a 'classicality' condition for many -body 
dynamics [45. ,46.1 . 

Here we extend the analytic discussion including a non- 
zero macroscopic field and apply it to a non-relativistic theory 
for an ultracold gas with, in the quantum case, bosonic statis- 
tics. For our comparison of quantum and classical dynamics 
we employ the functional-integral approach of the 2PI effec- 
tive action. We recall that the difference between the quantum 
and classical statistical theory can be expressed in terms of in- 
teraction vertex terms for the quantum theory which are absent 
in the classical statistical theory. As a consequence, the clas- 
sical generating functional is characterized by an important 
reparametrization property, which allows one to scale out the 
dependence of the dynamics on the scattering length a. As a 
consequence, for the classical dynamics the effects of a larger 
self-interaction can always be compensated by a smaller den- 
sity. It is shown that quantum corrections violate this invari- 
ance property. They become of increasing importance with 
growing scattering length or reduced density. This is used to 
derive a condition which may be used for experimenters to 
find signatures of quantum fluctuations when preparing and 
probing the dynamics of ultracold gases. This condition is not 
based on thermal equilibrium assumptions and holds also for 
far-from-equilibrium dynamics. As an application, we study 
the thermalization of a homogeneous one-dimensional ultra- 
cold Bose gas starting from a far-from-equilibrium initial state 
following Ref. [36]. We compare quantum and classical evo- 



lution and demonstrate the validity of our criterion for the 
nonequilibrium quantum field theory to be well approximated 
by its classical counterpart. To round off the analysis we dis- 
cuss some distinctive properties of quantum versus classical 
statistical dynamics. For instance, the decay of correlations 
with the initial state happens faster in the classical statistical 
theory for the one-dimensional Bose gas. Including quantum 
corrections the system remembers longer the details about the 
initial conditions. 

The use of functional methods to describe the dynamics of 
classical correlations dates back to the^ work of Hopf in the 
context of statistical hydrodynamics WA . A field theory for 
the description of classical fluctuations in terms of noncom- 
mutative classical fields was first suggested by Martin, Siggia, 
and Rose (MSR) ]48] and has been extensively used in critical 
dynamics near equilibrium [49il . This theory has been refor- 
mulated later in terms of Lagrangian field theory employing 
functional methods fs^ H M M 111 HI . In these field 
theoretical approaches to classical statistics, a doubling of the 
degrees of freedom occurs. For example, in the generating 
functional for Green's functions, besides each field appearing 
in the fundamental Lagrangian, a second 'response' field is 
integrated over The functional integral approach to quantum 
field dynamics developed by Schwinger and Keldysh employs 
a closed time path (CTP) contour ifSals? !] in the time-ordered 
exponential integral. The doubling of fields in the MSR and 
Lagrangian approaches to classical dynamics corresponds to 
the fields evaluated separately on the two branches of the 
Schwinger-Keldysh CTP lUllsllIM]- ImpHcations of the dif- 
ferences between the classical and quantum vertices, similar 
to the case considered in this article, have been discussed, for 
other theories, e.g. in Refs. S |6ll IH IH IH IH S . 

Our article is organized as follows: In Section |ll] we recall 
the functional description of quantum as well as classical sta- 
tistical non-equilibrium dynamics and use this to construct the 
respective 2PI effective actions. We then derive the time evo- 
lution equations and compare the nonperturbative expansion 
in the numbers of field components to next-to-leading order 
for the non-relativistic quantum and classical statistical the- 
ory. In Section ini] we present numerical results for the quan- 
tum and classical evolutions of an ultracold one-dimensional 
Bose gas. Our conclusions are drawn in Section |IV] In an 
appendix we provide some details of employed initial-state 
density matrices. 



II. CLASSICAL VERSUS QUANTUM DYNAMICS OF AN 
INTERACTING GAS 

We consider an ultracold gas of atoms with bosonic statis- 
tics. At sufficiently large phase-space densities the system can 
undergo a phase transition and form a Bose-Einstein conden- 
sate, given that dimensionality and trapping geometry fulfill 
the necessary conditions. For a dilute gas, i.e., if the atomic 
distance is much smaller than the characteristic length scale of 
the interactions, typically the s-wave scattering length a, the 
system may be described by a complex scalar field theory. 

We consider such a non-relativistic quantum field theory for 
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a complex-valued field (p{x). This fluctuating field is charac- 
terized by the Lagrangian density 



—diip*{x)di(p{x) - V {x)ip* {x)(p{x) 
Am 

i-{^*{xMx)r 



(1) 



in the defining functional integral for correlation functions as 
described below. We use units where h = 1. The space-time 
variable is a; = {xo = t,x.). It is summed over double in- 
dices i = 1, ...,d for d spatial dimensions and di = d/dxi. 
Here V denotes an external potential, and g a real-valued cou- 
pling constant. The Euler-Lagrange equation of motion de- 
rived from ([U reads 



id.^„ip{x) 



2m 



V{x)+gW{x)\^ 



^{x). (2) 



In the context of the physics of quantum gases of indistin- 
guishable Bosons, Eq. (|2]l is the well known Gross-Pitaevskii 
equation if the fluctuating field if{x)is identified with its 
quantum statistical average {f{x)) @ [l^. This equation 
approximately describes the time evolution of an inherently 
quantum system, a Bose-Einstein condensed ultracold gas, ne- 
glecting all quantum statistical fluctuations contributing to the 
dynamics of the scalar field. As such, the Gross-Pitaevskii 
equation is a classical field equation and the inclusion of quan- 
tum and statistical fluctuations beyond the classical field ap- 
proximation is described in the following. 

For this we recall in this section the field theoretical for- 
mulation of the many-body quantum and classical statistical 
time evolution. This gives that the dynamic equations for cor- 
relation functions differ only by certain terms in the quantum 
equations which are absent in the classical ones. As a side 
result, one recovers that the dynamics in the Hartree-Fock- 
Bogoliubov (HFB) rt67[|68ll69ll approximation is the same for 
quantum and classical statistical descriptions for same initial 
conditions. Moreover, the differences between the quantum 
and classical dynamical equations are identified in the non- 
perturbative 2PI 1/J\f approximation which goes far beyond 
HFB. 

While the field in Eq. ([T]) is assumed to be complex-valued, 
we will, in the following, switch to a representation of the field 
in terms of its real and imaginary part, ip = {ipi + iip2)/ V^, 
where the classical action reads 13611 
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(Pa{x)ipaix)(pbix)(pb{x), (3) 



with = J da;o / d'^a;. The free classical inverse propagator 
is given by 

^Go^abi^^y) = ^(2; - y) [~i<^lbdxo - HiB{x)Sab] , (4) 



where Hib{x) = —df /2m + V{x) denotes the single-particle 
Hamiltonian with interaction potential V{x) and 





i 



(5) 



the Pauli matrix in field-index space. Summation over double 
indices a,b,c. . . = 1, 2 is implied. Many of the following 
formal derivations are independent of the detailed form (HI 
of Gq^ and, in particular, equally valid for relativistic field 
theories. Note however that as Gq^ in Eq. ^ contains only a 
first-order time derivative, the canonically conjugate field, in 
the non-relativistic theory, is 



TTaix) 



6S[ip] 



S{dxa(Paix)) 



i<^lbVb{x), 



(6) 



in contrast to the relativistic case where the canonical momen- 
tum equals the time derivative of the field. 



A. Quantum statistical dynamics 

1. Generating functional 

For a given initial-state density matrix pnito), which may 
characterize a system also far from equilibrium, all informa- 
tion about the quantum field theory is contained in the gener- 
ating functional for correlation functions: 

Z[J,K;pd] = Tr[pD{to)Tcexp{t(^J^ J^{x)^a{x) 

^a{x)K^ab{x,y)My))} 



xy,C 



, (7) 



with Heisenberg field operators ^a{x) which obey for the 
non-relativistic theory the commutation relations 



[*a(i,x),$b(t,y) 



(8) 



In Eq. (I7]i, 7c denotes time-ordering along the closed time 
path C leading from the initial time to along the real time 
axis to some arbitrary time t and back to to, with ^ = 

dxQ J d'^x. Contour time ordering along this path corre- 
sponds to usual time ordering along the forward piece C+ and 
antitemporal ordering on the backward piece C^. Note that 
any time on is considered later than any time on C+. The 
source terms in Eq. O allow to generate correlation functions 
by functional differentiation such as 



(Tc$(xi)---$(x„)) = 



i"(5J(xi) • • •(5J(x„) 



(9) 



J,K=0 



where the field indices have been suppressed and we have 
used that for the closed time path Z = 1 in the absence of 
sources. We have introduced, in Eq. (|9]l, two contour source 
terms, J'^ and A''', which we use below to go over by Legen- 
dre transformation to the corresponding 2PI effective action. 



4 



2. Functional integral 

For a theory with action (O the generating functional 
Z[ J, K\ pj:,] can be expressed in terms of a functional integral 
using standard techniques (see e.g. Refs. ll44l [59tl and refer- 
ences therein): 



Z[J, K; pd] = J [d</7+] [d(^o ] po [(fi^ , 
X J Vip+Vif^ exp < z 

¥'+(*o,x)=ip,|(x) 



1 



xy 



a6 



(10) 



where poif^ ,(Po) — {(p~^\pD{to)\(p~) and the matrix el- 
ements are taken with respect to eigenstates of the Heisen- 
berg field operators at initial time, $a(ioi = 
(y9^^(x)|(y9^). In Appendix lAl we provide an explicit ex- 
pression for the initial-state density matrix poifQ , ^q) used 
later in our numerical calculations. The integral mea- 
sures are given as [dcp^] — Yia x '-^ft a(^) Vip^ = 
ria a;o>to x ^Vai^o, x), with the prime indicating that the in- 
tegration over the fields at initial time to is excluded. The 
superscript '+' ('— ') indicates that the sources are taken to 
be different on the forward (C+) and backward (C^) branch 
of the closed time path. Because of the different sources, the 
corresponding fields on the different branches are labelled ac- 
cordingly. The minus sign in front of the ' — ' terms accounts 
for the reversed time integration. Using this notation the ac- 
tion functional reads 



1 



iG-]., 



{flfl^tv't - fafafb'Pb) ' (11) 



which corresponds to the defining action (|3]l if the time in- 
tegration is replaced by an integration along the closed time 
contour C. 

In order to simplify the comparison with the classical sta- 
tistical field theory below, a standard linear transformation R 
of the fields is introduced as 



proliferation of symbols we have used here ipa, which agrees 
with the defining field in (O only for (p'^ — ip'^. Since this will 
be the case for expectation values in the absence of sources, 
where physical observables are obtained, and since there is no 
danger of confusion in the following we keep this notation. 
Correspondingly, we write for the source terms 



Ja 
Ja 



= R 



^ab ^ab 
^ab ^ab 



(14) 
(15) 



Inserting these definitions into the functional integral ( fTOb and 
using that ( fT4l ) and ( fTSl ) can be equivalently written as 



Ja 
Ja 



^ab ^ab 
^ab ^ab 



Jt 

-Ja 



-K, 



-K 



ab 
^ab 



(16) 



R-'^ (17) 



one finds: 

Z[jJ, K^,K^,K^,K^;pd] 

[d(^o][d<^o] PD [(Po + <^o/2, (po - (^o/S] 



J V'pV'>pexp\i S[p,^]+ J 



1 

xy 



{'Pa,<Pa) 



Ja 
Ja 



{Va,^a) 



(K^b 






{K^b 


K^bJ 





(18) 



where S[lp, 0\ = So[tp, ip] + Sint[p, 0] consists of the action 
for the free field theory 







So[^,^] = l I {^a,^a){J-, 
^ Jxy \ '^0,ab 



and the interaction part 



fpb 
Vb 



(19) 



Sint[f,0\ = ^1 / Va{x)(pa{x)pb{x)(pb{x) 



Va{x)Pa{x)pb{x)^b{x)- (20) 



3. Connected one- and two-point functions 



•Pa 



^a 



(12) 



From the generating functional for connected correlation 
functions 



where 



R = 



1 1 

2 2 
1 -1 



R- 



(13) 



such that Pa = {^t + t^a )/2 and (pa = pt -'Pa^ or, p+ = 
+ 'Pa/'^ and p~ ^ Pa — <^a/2, respectively. To avoid a 



W = -i InZ 
we define the macroscopic field (f>a, and (j>a by 



SJaix) 



6Ja{x) 



(21) 



SW 6W 

'l^o.ix), , =(j)a{x). (22) 
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The connected statistical correlation function Fab{x,y), the 
retarded/advanced propagators G^l^{x, y), and Fab{x^ y) are 
defined by 



5W 



dW 
^K^^,{x,y) 
5W 

SW 



i {(j)a{x)<i)b{y) +Fab{x,y)) 

\ {<l>a{x)My) - ^Gabi^,y) 

i (j>a{x)(f>b{y) - iG^h{x,y) 



- {(t)aix)(t)b{y) + Fab{x,y) 



(23) 



For vanishing sources the field = and the propagator 
F = lf59ll . Moreover, the retarded and advanced correla- 
tors G^{x,y) and G-'^{y,x) vanish for xq < yo- (See also 
the discussion in Sect. Ill B"3] ) In the absence of sources, these 
two-point functions are related then through the transforma- 
tion o, 



Fab 







'-'ab ^ab 



(24) 



to the time-ordered correlation functions written in the '±' 
basis. The inverse of the two-point function matrix (|24] | reads 







^iG^)ab 



X 



ab 



(25) 



where 

Kbi^^y) 



{G'^)-J{xrAFcd{z,w){G^)-^,\w,y) 
[(G«)-i.F.(GA)-i]^J:r,y). (26) 



The last equation introduces a compact matrix notation that 
will be employed below. 



4. 2PI effective action 

The 2PI effective action is obtained as the Legendre trans- 
form 

r = W~ f Ua{x)Jaix) + 4>a{x)Ja{x)) 

J X 

[K!b{x.y){ct>a{x)^b{y) + Fab{x,y)^ 

+ K^bix,y) (^Mx)My)~iG^b{x,y)) 

+ K^b{x,y) (^Mx)My)-iGtb{x,y)) 

+ K^bi^.y) (^Mx)My) + Fab{x,y)) }. (27) 

This corresponds to the 2PI effective action originally dis- 
cussed in Refs. ll3Tl[32i[33|] . however, written in terms of the 



(a) 



(b) 



x,a 



y,b 



FIG. 1 : (color online) (a) Diagrammatic representation of the corre- 
lators in the ip-if) basis. A full line indicates the 1- or lys-component, a 
broken line the 2- or (^-component. Fab{x, y) is the statistical corre- 
lation function, Gab{x,y) = pab{x,y)0{xo - yo) and Gab{x,y) = 
—pab[x, y)6{yo — xo) the retarded and advanced Green's functions, 
respectively. Their representation in terms of the real-valued spectral 
correlation function Pab{x, y) exposes the ^-functions which imply 
the respective time ordering in xq, yo. (b) Diagrammatic expansion 
of the quantum vertex term entering the action ([3}. The classical 
action i53\ is lacking the second contribution (red square). 



rotated variables. From ( |27l ) one observes the equations of 
motion for the fields 



sr 



S(t>a{x) 



-Ja{x)- / iK!,{x,y)4>b{y) 
+ \Ktb{^, y)My) + \My)Kt{y, x) ) , (28) 



Ha{x) 



-Ja{x)- j (Kf^,{x,y)Uy) 



+ ]^K^,{x,y)My) + \My)Kta{y.^)), (29) 



as well as for the two-point functions 



\Kj^b{x,y), 


(30) 


\Ktb{x,y), 


(31) 


\K^b{^,y), 


(32) 


\K^b{x,y)- 


(33) 



5Fab{x,y) 
5V 

''^Glix,y) 
ST 

''SG^x.y) 
6T 

SFab{x,y) 



The diagrammatic calculation of the 2PI effective action in- 
volves all closed two-particle irreducible Feynman graphs 
with fines associated to the full two-point correlators IBU [32l 
I33I1 . This is illustrated in Fig. [T] for vanishing macroscopic 
field (j> and in the absence of external sources. 

In the presence of a non-vanishing field value, 7^ 0, 
the interaction vertices are obtained from (|20l) by shifting in 
S'[(,c, (^] the field (p ^ (j> + (p, and collecting all cubic and 



quartic terms in the fluctuating fields and (p, i.e. 

5'int[</3, <^;<?i'] = -| J'fa{x)fa{x)ipb{x)ifb{x) 
-| j (pa{x)ipa{x)(pb{x)ipi,{x) 

J X 

[ 'fa{x)<j)a{x)fb{x)iph{x) 

>fa{x)^Pa{x)>fb{x)<j)b{x). (34) 



The quadratic terms in the fluctuating fields are taken into ac- 
count in the classical inverse propagator (01) by the replace- 
ment 

HiBSab —> HiB + ^<j)c{x)<j)c{x) Sab + ff'/'a (a;)0b (x) . (35) 

In the free part of the action (T% . and in the dynamic equa- 
tions derived below, this corresponds to a field dependent 
iGQ^{x,y;(l)), while the general form of the equations re- 



mains unchanged 117411 . We note that linear terms in the fluctu- 
ating fields ensure cancellation of possible tadpole contribu- 
tions, which therefore do need not to be considered explicitly. 



5. Exact evolution equations 

For vanishing sources, the exact inverse two-point function 
can then be written as l31. ,3Z,,33il 

z{G^): ' 

KG^rab Kb 

G, 
Ga.lb 



Q,ab 





-^^tb 
yF 
^ab 



(36) 



where X is defined in Eq. (|26] | and the retarded, advanced 
and statistical self-energies are related through the transfor- 
mation ( fT3T l. 



= {R-'f ( -E"^ ) R-\ 

\ ^ab ^ab J 

(37) 

to the self-energies written in the '±' basis. To the re- 
tarded/advanced self-energy E^/^ and the statistical self- 
energy contribute only graphs with propagator lines as- 
sociated to G^'^ and F, which can be obtained from closed 
two-particle irreducible graphs by opening one propagator 
line iMmill. 

To convert Eq. (l36T l for the inverse propagator into an equa- 
tion, which is more suitable for initial value problems, we 
convolute with the propagator matrix (|24] | from the right and 
with the classical propagator from the left. This yields the 
Schwinger-Dyson equations for the retarded/advanced propa- 
gator in the absence of sources. 



^R/A _ (^R/A _ (^R/A ^ j]R/A (^R/A 



(38) 



and the statistical propagator, 

F = FQ^Fn-Y.^ -G^ -G^ ■ [E^ -F + yF ■G^\, (39) 

where we have used the compact notation introduced in 
Eq. ( |26] l. For the spectral function. 



Pab{x,y) = G\(x,y) - G^i,{x,y), 



(40) 



the Schwinger-Dyson equation follows from Eqs. ( |40l i and 

as 



P 



Po-G^ -Y^ -G^ +G^ -Y^ -G^. (41) 



Acting on Eqs. ( |39] l and (ITTl i with Gq^ from the left brings 
these equations in a form which is more suitable for initial- 
value problems. For this we write 

Sab(a;,y;0) = ^fb ix)Six - y) + e{xo - yQ)T,P^^{x,y; (f)), 

(42) 

where we introduce the spectral part of the self-energy 

y; 0) = J:^,{x, y; 0) - y; 0). (43) 

Using that Fq and po satisfy the homogeneous field equations 



Go,Li^,z)FQ^cb{z,y) 



0, 



(44) 



and taking into account all 0-functions one obtains the dy- 
namic equations for the two-point correlation functions F and 

pM^m-- 

[i'^Ldxo + Mac{x)] Fcb[x, y) 



dz j:P^{x,z:<j>)Fcb{z,y) 



<izT,^cix,z;4))pcbiz,y), 



[if^lAo + Mac{x)] Pcb{x,y) 



dz J:P^^{x,z;(j))pcbiz,y), 



(45) 



where we employ the notation 



/ dz= dzo / d'^z. (46) 

J to tQ J 



Here 

Mab{x) 



= Sat 



H 



IB 



{x) + I [(i)c{x)(j)c{x) + Fcc{x, X)^ 



+g(^(j)aix)4>b{x) +Fabix,x)^ 



(47) 



is the mean-field energy matrix which includes the field- 
dependent terms of the classical inverse propagator ^Gq 
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defined in Eqs. ^ with (1351 ), and the local part of the self- 
energy, 

^abi^) = ^SabF,,{x,x)+gFabix,x), (48) 

as defined in Eq. ( |43] l. We note that the spectral function at 
initial time is characterized by the commutator (O, with 



(49) 



for the non-relativistic theory. 

The dynamic equation for the mean field (j)a (x) is obtained 
from Eq. ( |29] l with (i>a{x) = in the absence of sources. It 
reads 1 35,. ,36,1 : 



( - i'^lb^xo - gFab{x,x)^(l)b{x) - (^Hib{x) 

+ ^[(l)cix)(t)cix) + Fccix,x)A4>aix) = 



5(j}a{x) 



(50) 



where the functional derivative is taken for fixed two-point 
functions and r2 contains all closed 2PI graphs |„3 1 . 32, 33J- 
In Sec. Ill CI we will employ an expansion in the number of 
field components to next-to-leading order to approximately 
describe the dynamics, for which one obtains the compact ex- 
pression issli 



5Vn 



S<j)a{x) 



dyi:P,{x,y;cl>=0)My)- (51) 



B. Classical statistical dynamics 

1. Classical equation of motion 

The classical field equation of motion can be derived from 
the defining action Q. In the basis (fT2l i the same equation 
of motion for the field ipa (x) is obtained from the action with 
dull and (m by 



dS[^, If] 



5(p 



0, (52) 



as a consequence of setting = 0. For the first identity in ( 
we have defined 5''^' [</?, (p\ = [</3, (^] + S^l^ [(^, (p\ as that part 
of SYp, (p\ which is linear in (p. As a consequence, the non- 
interacting part 5*0 Yp, ip] agrees with the respective expression 
appearing for the quantum theory (fT9l) . The interaction part of 
S"^'[(p, 0\, however, differs from ( |20] | in that it contains fewer 
vertices: 



SttVf, "^l = ~f / ^a{x)(Pa{x)ipb{x)ipb{x). 



(53) 



This is illustrated in Fig. [TJ). The absence of vertices beyond 
those which are linear in (p turns out to be the crucial differ- 
ence between a classical and a quantum statistical field theory 
as is discussed in the following. 



2. Classical statistical generating functional 

We will construct the generating functional for the classi- 
cal statistical field theory, Z'^\j, K; W] for given probability 
functional W for the fields at initial time, similar to the ex- 
pression (fTsT l for the quantum theory. For this, we rewrite the 
equation of motion ( |52] i as a (5-constraint in a functional inte- 
gral using the Fourier transform representation 



5(p 



^ J Dip exp^^iS''\ip,'p]+i ^ 7ro,a(x)<^o,a(x)| . (54) 

For the last equality, the second term in the exponent subtracts 
the boundary term at the initial time to, which follows from 
partial integration. We emphasize that we have to take into ac- 
count this boundary term since to will be taken to be finite for 
the nonequilibrium initial-value problems considered. Note 
that in the path integral (fTsT l, the action S[(p, 0\ is implied not 
to depend on boundary terms at initial time Iq. (As a conse- 
quence, the term i ttq ,j(x)(^o,a(x) appears in Eq. ( |54l i and 
not half of it.) Here 7ro,a(x) = T^aixo = to,yi) is the initial 
canonical momentum, with TTa{x) = ia'^i^iptix) for the non- 
relativistic theory as in Eq. (|6]l. We note that Eq. (l54l i remains 
valid also for a relativistic theory, where 7ro.a(x) denotes the 
time derivative of the field at time to. 

It should be stressed that the Fourier transform expression 
( |54l i could not be achieved with S'^^[(p, ip] replaced by S[(p, ip] 
since the latter is not linear in ip. For the non-interacting the- 
ory, however, Sf^[(p, ip] = So[(p, ip] holds, and we will recover 
the fact that the free classical and quantum theories are gov- 
erned by the same dynamics for same initial conditions. The 
same fact holds true for all Gaussian (leading-order large-A/^ 
or Hartree-type) approximations. Therefore, it is crucial for a 
quantum classical comparison that we go beyond and consider 
the next-to-leading order corrections as we do in Sect. Ill A 41 
below. 

Any functional f[f'^^] of the classical field ip"^^, being a so- 
lution of the classical equation of motion ( |52] i, can be written 



-SS<='[^,ip] 



V'^f[ip]5 



5(p 



X exp 
where the Jacobian reads 



cl 


\lS''\ip,(p\+l j 7ro,a(^0,a}j^M, (55) 



det 



5^S'^\ip,^] 
S(pS(p 



(56) 
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Here the Jacobian plays the role of an irrelevant normalization 
constant, which has been discussed in detail in Ref. [66] and 
references therein. 

Classical correlation functions are obtained as phase-space 
averages over trajectories given by solutions of the classical 
field equation (|52|. Such averages, for an arbitrary functional 
of the field, are defined as 



(57) 



Here W^[93q', 7rg'] denotes the normalized phase-space proba- 
bility functional at initial time (and is not to be confused with 
the generating functional (ISTTl). In Appendix|A]we provide an 
explicit expression for T/F[i^q', ttq'], a functional of four fields 
which is symmetric under the exchange of ttq'^ and i<Jab'p'ob 
for a = 1,2, since the canonical momentum (|6]l is propor- 
tional to the field itself. The measure indicates integration 
over the classical phase-space. The theory may be defined on 
a spatial lattice to regulate the Rayleigh- Jeans divergence of 
classical statistical field theory. 

Using Eqs. ( |55T l and (|57T i. we can now write down a gener- 
ating functional for classical correlation functions in the form: 

Z^V, J, K^, K^,K^; W]= [ [d^S'][d^§']W^[^S\ 



/ T)' tpDip exp < i 



Ja 
Ja 



^ab ^ab 
^ab ^ab 



ffb 



(58) 



where S"^'[(/9, (^] = So[(p, (f]+Sf^^^[(p, (f] as defined in Eqs. (fT9]l 
and ( l53T l, and where the relation (|6]) between the field and the 
canonical momentum is implied. 

Comparing with the quantum generating functional in 
Eq. (fTSl l. and using that the Jacobian J^[(p] plays the role of 
an irrelevant normalization constant, we find that the classical 
functional ( fSSl l takes the same Lagrangian form, with initial 
conditions given by a density matrix which is character- 
ized by the Fourier transform of the phase-space probability 
distribution 11^ [(/Sq, ttq]: 



Pd [vq + <^o/2, V3o - "^0/2] 

= / [d7ro]T4^[v3o, ti'o] exp < i / ttq^q (^o,a 



(59) 



In summary, the generating functionals for correlation func- 
tions are very similar in the quantum and the classical statis- 
tical theory. The crucial difference is that the quantum theory 
is characterized by more vertices. As a consequence, we can 
follow the very same steps as in Sect. Ill Al to construct the 
classical statistical 2PI effective action and to derive the time 
evolution equations from it. 



Similar to the discussion for the quantum theory above, for 
(j> ^ the interaction vertices for the classical statistical the- 
ory are obtained by shifting in S^^[(p, if] the field (p (p + ip, 
and collecting all cubic and quartic terms in the fluctuating 
fields, i.e. 



g / (paix)ipa{x)(pb{x)(l)bix) 



- I Lpa{x)(t)a{x)Lpb{x)ipb{x). (60) 



This can be compared to the respective expression for the 
quantum theory, Eq. (l34l i. 



3. Correlation functions 

Comparing quantum and classical statistical dynamics 
amounts to comparing the time evolution of classical statisti- 
cal n-point functions with that of the respective quantum ones. 
The classical functions are obtained as phase-space averages 
(|57] | over trajectories given by solutions of the classical field 
equation (l52b . As an example, the macroscopic or average 
classical field is given by 



(lf^[x) = (^a(x))e 



(61) 

This is equivalent to the field average obtained from the gener- 
ating functional ( |58l l by functional differentiation with respect 
to Ja{x). Similarly, also (t>'^a{x) can be obtained by functional 
differentiation with respect to J a [x) in the same way as de- 
scribed in Sect. Ill A 31 if the quantum generating functional 
in Eq. ( 1211 1 is replaced by its classical counterpart defined in 
Eq. (ISSl l. Moreover, as for the quantum case, 4)a{x) = for 
vanishing sources. This can be directly seen by taking the 
functional derivative of Z^^ with respect to Ja{x), setting all 
sources to zero except Ja{x), and performing all steps back- 
ward in Eqs. dSSl l and (l55T l. This is possible as the remain- 
ing source term is linear in ip. Using the same reasoning one 
finds that also the two-point correlation function Fab{x,y), 
which is defined as in Eq. ( |23] l as a functional derivative of 
W^' = —i InZ'^' with respect to K^, vanishes identically in 
the absence of sources. 

Similarly, defining the connected classical statistical prop- 
agator Fl^l{x, y) according to Eq. ( |23] l results in 



(62) 



The classical equivalent of the quantum spectral function 
Pabi^> y) obtained as follows: Taking the functional deriva- 
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lives according to Eq. ( |23] ) leads, with Eq. ( |40] |. to 



X I 



X exp ■ 



S'^^ifi'f] + / T^Q.a'fQ.a 



(63) 



If Xq < yo the expectation value of Lpa{x)(pb{y) vanishes fol- 
lowing the same reasoning as for the field (jf^, in accordance 
with the fact that this term corresponds to the retarded propa- 
gator G^(a;, y). Analogously, if xq > yo, the advanced prop- 
agator, i.e. the average of (pa{x)ifb{y), vanishes. 

We consider the case that x^ — tg. If yo > xo, one obtains 
from Eq. ( l63T l and the definition (|6]l of the canonical field mo- 
mentum: 



X exp < z 



[d^^oWKW[^o\< 



1 /• Va'(io,x) 5^f{y) 



(64) 



Extending this procedure to to ~ yo < xo one recovers 
the well-known fact that the classical spectral function is ob- 
tained by replacing — i times the commutator with the Poisson 
bracket: 



pfb{x,y) ^ -({(/3a(a;), Vfc(y)}pB)ci- 

The Poisson bracket with respect to the initial fields is 



(65) 



{(Paix),ipbiy)} 



PB 



(57ro,c(z) Sipo.ci^) 



(66) 



(Summation over c). Note that, in order to arrive at Eq. ( |65] ). 
we used that the Poisson brackets are invariant under the clas- 
sical time evolution of the fields and therefore valid for any 
times Xo, yo- 

As a consequence, one finds the equal-time relations for the 
classical spectral function: p'^i,i^,y)\xo=yo — ~^'^ab^i^ ~ 
y). Note that they are in complete correspondence with 
the respective quantum relations in Eq. (|49] l. Equivalently, 
the free spectral function Po'cbl^'?/) '^^'^ statistical function 
-^0 cb('^' y) solutions of the homogeneous equations cor- 
responding to Eq. (l44b . with initial conditions determined for 
Pq' by the equal-time canonical relations, and for Fq' by the 



initial probability functional W^[(/3oj '''o]- Also for the classi- 
cal statistical theory Eqs. ( |50] l and ( |45] l are the exact equa- 
tions for the field cj) and the correlation functions F and p, 
respectively. There is a difference between the classical and 
quantum equations of motion only in the self-energy contribu- 
tions corresponding to and S''. This difference arises from 
the different properties of the interaction part of the quantum 
and classical actions ( |20] | and ( 153] ), respectively. We discuss 
this difference for the non-perturbative 2PI 1 /Af expansion to 
next-to-leading order in the following subsection. 

Summarizing, one finds when comparing classical statis- 
tical and quantum many-body dynamics that the generating 
functionals for correlation functions are very similar in the 
classical and the quantum theory. However, the quantum the- 
ory is characterized by more vertices. As a consequence, the 
same techniques can be used to derive time evolution equa- 
tions of classical correlation functions that are employed in 
quantum field theory. Eventually, in the basis corresponding 
to the fields (p and if, the classical dynamic equations have the 
same form as their quantum analogues but are lacking certain 
terms due to the reduced number of vertices. 



C. Quantum versus classical statistical evolution 



The classical statistical generating functional ( 1581 ) exhibits 
an important reparametrization property: If the fluctuating 
fields are rescaled according to 



(67) 



then the coupling g drops out of 5'^'[(y5, i^] = So[ip,0] + 
SfjJt [(fi, (p] defined in Eqs. (fT9] l and (53[ . The free part So [ip, p] 
remains unchanged and the interaction part becomes 



p',{x)ip',{x)p',{x)p',{x). (68) 



Moreover, the functional measure in Eq. iSEi is invariant un- 
der the rescaling ( |67] |. and the sources can be redefined ac- 
cordingly. Therefore, the classical statistical generating func- 
tional becomes independent of g, except for the coupling de- 
pendence entering the probability distribution fixing the ini- 
tial conditions. Accordingly, the coupling does not enter the 
classical dynamic equations for correlation functions. All the 
g-dependence enters the initial conditions which are required 
to solve the dynamic equations. 

In contrast to the classical case, this reparametrization prop- 
erty is absent for the quantum theory: After the rescaling (|67| | 
one is left with S[(p', p'] whose coupling dependence is given 
by the interaction part 



'5'int [V'^'P'] 



'P'aix)v'aix)'P'bix)'P'bix) 



^ / ip',{x)^',ix)ip'bix)p'i,ix), (69) 
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according to Eq. ( |20] |. Comparing to ( |68] ) one observes that the 
'quantum' vertex, which is absent in the classical statistical 
theory, encodes all the g-dependence of the dynamics. 

The comparison of quantum versus classical dynamics be- 
comes particularly transparent using the above rescaling. The 
rescaled macroscopic field and statistical correlation function 
are given by 

'l^'ai.x) = Vg<l>a{x) , F^^{x,y) = gFab{x,y) , (70) 

while the spectral function Pab{x; y) remains unchanged ac- 
cording to Eqs. ( |23] | and ( l40b . Similarly, we define for the 
statistical self-energy entering the dynamic equations i 



1. Quantum versus classical statistical self-energy 



(71) 



To identify the precise difference between the quantum and 
the classical time evolution, details about the self-energies are 
required. In the following we will ernploy the 2P1 1 /M ex- 
pansion to next-to-leading order ll34l [ssll . This is a nonper- 
turbative expansion in powers of the inverse number of field 
components M which, in the context of a non-relativistic Bose 



gas, \& Af = 2 as discussed in detail in Ref. 
the result for the self-energies for Af = 2 ll34ll 35 
the rescaled variables (ITOb and (TtTI ) read: 



. We quote 
1, which for 



-[l'p{x,y)<l>'aM{y) 
+ [l'Ax,y) + P'F{x,y)]F:,,ix,y) 

2 

[Ipix,y) + Ppix,y)] pabix,y)y 

Kbi^^y) = -[lpix,y)(l}'aix)(l)'hiy) 

+ [Ip{x,y) + Ppix,y)]F^,{x,y) 
+ [lFix,y) + Ppix^y)] Pabix,y)y (72) 
The functions fp and Ip satisfy 

I'pi^^y) = -^'pi^^y) 

+ / dzlp{x,z)ll'p{z,y) 
Jta 



ipi^^y) 



with 



n'Fi^,y) 



rvo 

- / dzFp{x,z)np{z,y), 
-^pi.x,y) 

+ / dzlp{x,z)'\lp{z,y), 
•I yo 



FUx,y)FU{x,y) 



-^Pab{x,y)pabix,y) 



(73) 




(C) 



(Ql) 



(Q2) 



FIG. 2: (color online) Diagrammatic representation of the bubble 
chains contributing to the functions I'p and Ip at next-to-leading or- 
der in the 2PI 1 /f\f expansion. The meaning of lines and vertices is 
explained in Fig.[T] The chains of type (Ql) and (Q2) only appear 
for a quantum system. Type (C) is present both in quantum and clas- 
sical systems. These diagrams exhibit that in each term contributing 
to the functions I'p and Ip, there is at most one loop involving two 
correlators F' or p. In the classical limit of quantum theory, the loops 
involving two spectral functions p are suppressed compared to those 
with two statistical functions F' . 



The functions Pp and Pp, which vanish if cj) — 0, will be 
discussed below. Iterating the integral equation ( l73T l implies 
that the functions I'p and Ip can be represented as a series of 
bubble chains, as shown in Fig. |2] A full line represents a 
statistical correlator, while a retarded (advanced) propagator 
line changes, from left to right, from full (broken) to broken 
(full), see Fig. [U 

The quantum and classical vertices are depicted analo- 
gously in Fig. [Tfb). Connecting the correlators through the 
respective vertices, one finds which types of bubble chains 
appear. The classes of non-vanishing bubble chains shown 
in Fig. |2] confirm the structure of the functions I'p{x, y) and 
Ip{x, y) which are determined by the integral equation ( |73] i: 
In each term contributing to the diagrammatic expansion of 
these functions there is at most one loop containing two F' 
or two p correlators, the latter resulting from either or 
G^. In addition to this one finds that only the loop contain- 
ing two p correlators goes with the vertex which is present 
in the quantum case only. Hence, considering the classical 
dynamics, the ^ g^p^ terms are absent together with the 
'quantum' vertex. For the same reasons, the contribution 
~ (5V4)[^p(a;, y) + Pp{x, y)]pab{x, y) to T,^f^{x, y) is absent 
for the classical statistical theory. 

We proceed by considering the functions P'p and Pp, which 
are relevant for the case of a non-vanishing 0'. The classes of 
bubble chains appearing in the self-energy for 7^ are de- 
rived from those shown in Fig.|2l by replacing at most one full 
F^ij{x,y) line by the product of mean fields. This 

means, that the pair of correlators in one loop are replaced by 
either of the combinations 

H'pix,y) = -4>'a{x)F'^b{x,y)4>'b{y), 

Hp{x,y) = -(j)'a{x)pab{x,y)(l)'b{y). (75) 



^pi.x,y) = -F^i,{x,y)pab{x,y). 



(74) The functions P'p and Pp entering the self-energies ( f72] i then 
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read S: 



— / dz 


\Hp{x, z}Ip{z, y) + Ip{x, z)Hp{z, y}\ 


J to 




rvo 




+ / dz 


[Hp{x,z)Ip{z,y) + rp{x,z)Hp{z,y)] 


J to 




rxo 


rVo 


- dz 


I dvIp{x,z)H'p{z,v)Ip{v,y) 


J to 


Jto 


rxo 




+ / dz 


/ dv Ip{x,z)Hp{z,v)I'p{v,y) 


J to 




rVo 




+ dz 


[ dvI'p{x,z)Hp{z,v)Ip{v,y)y 


J to 


J Zo J 


Pp{x,y) = - 


-gl^Hp[x,y) 



rxo 

- dz [Hp{x,z)Ip{z,y) + Ip{x,z)Hp{z,y)] 
Jyo 

rxo rzo ~| 

+ / dz / dvIp{x,z)Hp{z,v)Ip{v,y)\. (76) 



We observe that there are no differences in the functions P'p 
and Pp in the classical statistical limit except for a dependence 
on modified I'p and Ip, since the terms involving the mean 
field (f)'^ correspond either to [F'Y or to F' p loops. This im- 
plies that the presence of a non-vanishing mean field 0' does 
not add characteristic 'quantum' terms to the dynamical equa- 
tions of motion. 

Moreover, the tadpole diagrams contain the classical vertex 
and the F' correlator only. This shows that the left hand sides 
of the quantum dynamic equations, Eqs. (ISOl l and (|45] |. do not 
contain any contribution proportional to the quantum vertex 
marked with a (red) square in Fig. [Ub). As a result, the dy- 
namics in the Hartree-Fock-Bogoliubov ll67l IHsi l69ll approxi- 
mation is the same for the quantum and the classical statistical 
theory for same initial conditions. Differences arise, at most, 
in the self-energies S^' and Yi^ . We point out that, neglecting 
the right-hand sides of Eqs. (ISOl l. (l45l l. these equations consti- 
tute a set of time-dependent HFB equations for the mean field 
and the two-point functions, cf., e.g. Refs. lisjlfSTll . In this ap- 
proximation, the dynamics of p decouples from that of cj) and 
F. Neglecting also F, Eq. ( l50b becomes the Gross-Pitaevskii 
equation. Quantum fluctuations play no role in these approx- 
imations. For this reason, HFB is commonly termed a mean- 
field theory ill . 

Summarizing, one concludes that all equations ( |45] l, ( ISOl l, 
( |47] i, ( |72] ), ( |73] l, ( |74l ), ( |75] l, and remain the same in the 
classical statistical limit except for differing expressions for 



the statistical components of the self-energy 

yF'iT ,,^ classical limit / , , CrU' fl/^ 

+ [I'p{x,y) + P'p{x,y)]FU{x,y)^, 

U'p{x,y) ^lF:^,ix,y)FU{x,y), (77) 

replacing the respective expression in Eqs. ( |72] ) and ( l74l i. The 
classical statistical self-energies can be obtained from the re- 
spective quantum ones by dropping two spectral (p-type) com- 
ponents compared to two statistical (F-type) functions. For 
vanishing macroscopic field^ where Pp^ = this corre- 
sponds to the result of Ref. |45 ll . Using the rescaled variables 
( ITOI ) and ( fTTT i the 'quantum' terms can be directly identified 
since they are the only g-dependent terms, which are absent in 
the classical statistical theory according to the above discus- 
sion. As a consequence, for the classical dynamics the effects 
of a larger coupling can always be compensated by changing 
the initial conditions such that Fg, as well as 4>y/g, remain 
constant. This cannot be achieved once quantum corrections 
are taken into account, since they become of increasing im- 
portance with growing coupling or reduced initial values for 
F and 0. 

2. 'Classicality' condition 

Eq. ( ITTI i describes the differences between quantum and 
classical statistical equations of motion. In turn one can ask 
under which conditions these differences are negligible. In 
that case the quantum dynamics can be well approximated by 
classical statistical dynamics. To analyse this we iteratively 
expand Eq. ( |73] ) in terms of li'p and lip and compare term by 
term the statistical components of the quantum self-energies 
in Eqs. ( |72] | and ( |74] | to the respective classical ones according 
to ( ITT] ). One finds that a sufficient condition for the suppres- 
sion of quantum fluctuations compared to classical statistical 
fluctuations is given by 

3 

\Kb{x^y)Fcd{z,w)\ > -5^ \pab{x,y)pcd{z,w)\ . (78) 

This condition is not based on thermal equilibrium assump- 
tions and holds also for far-from-equilibrium dynamics. In 
particular, it is independent of the value of the macroscopic 
field (f). Our results, therefore, agree with previous i/i = es- 
timates in Refs. [44, 45.1 ■ We also note that the condition ( |78] l 
is precisely the same as the one obtained from the 2P1 loop 
expansion fsi','^?]. 

The condition ( |78] l can be applied, of course, also in ther- 
mal equilibrium, for which the statistical and spectral corre- 
lation functions are no longer independent quantities. In ther- 
mal equilibrium they are rather related to each other through 
the fluctuation-dissipation relation, which, for a homogeneous 
system in energy-momentum space, reads 

p) = Q + n{uj, T)\ p^^^^ [oj, p), (79) 
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with the Bose-Einstein distribution function n{uj, T) = 
^^{uj-p,)/kBT _ For large temperatures, /cbT' lu — fi, 
one has |F('='i)'(w, p)|/5 > |p('''i)(a;, p)|, i.e., the classicality 
condition is fulfilled for all modes whose occupation number 
^ p)/g is much larger than 0(1). The equivalent 

statement can be directly derived from (fTSl l also for nonequi- 
librium evolutions whenever it is possible to define a suitable 
'occupation number' from a space-time or energy-momentum 
dependent proportionality between F and p. A nonequilib- 
rium example will be discussed in Sect. |III] 

Away from equilibrium the situation is often considerably 
more complicated. Strictly speaking the condition ( fTsT l must 
be valid at all times and for all space points, or momenta in 
Fourier space, for the classical and the quantum evolution to 
agree. In practice, however, it needs only be fulfilled for time 
and space averages. In Sect. |III] we will demonstrate how 
quantum evolution can be approximated for not too late times 
by classical statistical dynamics, if the correlation functions 
satisfy (iTST i at initial time. In order to have quantum fluctu- 
ations playing a significant role, also for dynamically evolv- 
ing gases, one either needs to increase the interaction strength 
g accordingly or change the phase-space structure by chang- 
ing the external trapping potential. For example, in a one- 
dimensional trap, an effectively strong coupling and strong 
quantum fluctuations can be induced by reducing the line den- 
sity of atoms while their interaction strength is kept constant. 
Such a case will be considered in Sect. 1111 Cl below. 



III. FAR-FROM-EQUILIBRIUM TIME EVOLUTION OF 
AN ULTRACOLD BOSE GAS 

In this section we apply the theoretical methods summa- 
rized above to describe the equilibration dynamics of a uni- 
form ultracold gas of bosonic sodium atoms which are con- 
fined such that they can move in one spatial dimension only. 
With present-day experimental technology, such a situation is 
achievable, e.g., using strong transversal confinement in an 
optical lattice or in a microtrap on the surface of a chip. 

We will compare the evolution involving only classical sta- 
tistical fluctuations with that which also takes into account 
quantum corrections. Both, the classical and the quantum 
gas are assumed to be initially characterized by the same ini- 
tial conditions far from thermal equilibrium. For the quantum 
gas, the ensuing equilibration process is found to happen on 
two different time scales. A fast dephasing period leads to 
a quasistationary state which shows certain near-equilibrium 
characteristics but is still far from being thermal. After this, 
the system approaches, within an at least ten times longer pe- 
riod, the actual equilibrium state. On the contrary, the clas- 
sical gas does not show the dephasing when considering the 
same initial mode occupation numbers as for the quantum gas. 
We show that the dephasing in the short-time evolution of the 
quantum gas can, to a certain extent, be simulated by a classi- 
cal gas if one chooses the same initial values for the correla- 
tion functions. Our results show explicitly that quantum fluc- 
tuations only play a role for modes whose occupation number 
is sufficiently small. Hence, the dynamics of the weakly in- 



teracting one-dimensional gas is almost purely classical. The 
long-time evolution is different in the classical and quantum 
cases, since only the latter can reach a Bose-Einstein distri- 
bution. We then study the example of a strongly interacting 
one-dimensional gas. For such a gas we find characteristic 
differences. For instance, for given identical initial values for 
the correlation functions, in the quantum evolution the decay 
of correlations with the initial state takes much longer as it 
would be expected from a calculation in the classical approx- 
imation. 



A. Initial conditions 

The 2PI effective action approach is convenient for situ- 
ations, where at time t ^ one has a Gaussian state, i.e., 
a state, for which all but the correlation functions of order 
one and two vanish 17511 . In the following we will consider 
a one-dimensional uniform system, for which the two-point 
functions Fab{x, y) and Pab{x, y) are spatially translation in- 
variant. We will therefore work in momentum space, where 
the kinetic energy operator is diagonal. Moreover, we choose 
the mean field (j) to vanish initially. Then, for reasons of num- 
ber conservation, the equations of motion ( fSOl ) and ( l45b will 
conserve (j) — for all times. We note that, since there is no 
spontaneous symmetry breaking in one spatial dimension at 
non-zero temperature, the field always approaches zero even- 
tually, irrespective of its initial value. 

Having prescribed initial values Fab(0,0;p), with 
Pafc(0,0;p) given by Eq. ( |49] ), the coupled system of 
integro-differential equations ( |45] ) yields the time evolution 
of the two-point functions, in particular, of the momentum 
distribution 

n(t,p) = ^(^Fn{t,t;p) + F22{t,t;p) - ay (80) 

For the quantum gas, one has a — 1 from the Bose commu- 
tation relations, while, for a gas following classical statistical 
evolution, a = Q. 

We choose, at f = 0, a Gaussian momentum distribution 

n(0,p) = ^e-^''/'^^ (81) 

which constitutes a far-from-equilibrium state for the quantum 
gas as well as for the classical gas if the interactions are non- 
zero and the corresponding interaction energy is much larger 
than the kinetic energy. 

The initial pair correlation function vanishes, 

0= ^(^Fii{t,t-,p) ~ F22{t,t;p)^ +iFi2{t,t-,p), (82) 

for t = 0,m accordance with total atom number conservation 
at non-relativistic energies |76]. Hence, 

Fn{0, 0; p) = F22 (0, 0; p) = n{0,p) + a/2, (83) 
Fi2(0,0;p) = F2i(0,0;p) = 0. (84) 
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As far as the spectral functions are concerned, in the quan- 
tum case, the Bose commutation relations, and, in the classical 
case, the Poisson brackets require, cf. Eqs. dUll, dH: 

pn{t,t;p)^P22{t,t;p) = 0, (85) 
-pi2it,t;p) ^ P2iit,t;p) = l. (86) 

We have investigated the dynamic evolution of a ID Bose 
gas of sodium atoms in a box of length L — NsUg, with pe- 
riodic boundary conditions. We choose the numerical grid 
such that it corresponds to a lattice of Ng points in coordinate 
space with grid constant as, and the momenta on the Fourier 
transformed grid are p„ — (2/as) sin(n7r/iVs). The results 
presented in the following are obtained using Ng = 64 modes 
on a spatial grid with grid constant = 1.33 ^m. We first 
consider a line density of the atoms in the box of rii = 10^ 
atoms / m. In this case the atoms are weakly interacting with 
each other, such that (7iD = ^1^7711/771, with the dimensionless 
parameter 7 — 1.5 • 10~^. The width of the initial momen- 
tum distribution is chosen to be ct = 1.3 • 10^ m^^. In order 
to explore a strongly interacting gas we then reduce the total 
number of atoms by a factor of 100, increasing the dimen- 
sionless interaction parameter 7 by lO''. Here it is important 
that we employ a nonperturbative approximation which is not 
based on weak interactions. 
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FIG. 3: (Color online) Momentum-mode distribution n{t; p) for the 
initial state (red filled circles, interpolated by red solid line) and 6 
subsequent times t until no change can be observed for t > 0.6 s. 
The interpolation of the final distribution is shown as a black solid 
curve. Note the logarithmic scale. The occupation number are nor- 
malized by the total number of atoms in the box, niL = 853. The 
gas is in a far-from-equilibrium state initially, characterized by a 
Gaussian distribution n(0; p), Eq. J81b. with width a = 1.3-lO^m^^. 
It is weakly interacting, 7 = 1.5 • 10"'^. Since we consider a homo- 
geneous gas and a symmetric initial state, the occupation numbers 
are invariant under p —p. 



B. Equilibration of the quantum gas 

To solve Eqs. ( |45] ), with the self-energies given by 
Eqs. ( |72] i, for the initial conditions given in the previous sec- 
tion, we have implemented a parallelized Runge-Kutta solver 
and used a cluster of 3 GHz dual processor PCs with up to 
one node per momentum mode. The correlation functions 
Fab{t,t';p) and pab{t,t';p) were propagated, for fixed t', 
along t, using a second-order Runge-Kutta algorithm. After 
each Runge-Kutta step, the Ip^p integrals were updated ac- 
cording to Eqs. ( |73T l. The dynamic equations derived from 
the 2P1 effective action are, by construction, number and en- 
ergy conserving. While number conservation, by virtue of the 
0(2) symmetry of each diagram, is given exactly, energy con- 
servation may be violated by the chosen discretization along 
the time axis. Hence, in order to ensure optimal energy con- 
servation numerically, a fourth-order Runge-Kutta algorithm 
was employed for the propagation of the correlation functions 
along the diagonal t ^ t'. 

Fig. [3] shows, as a (red) solid curve, the initial Gaussian 
momentum distribution of the gas, on a logarithmic scale, 
where it forms an inverted parabola. The filled circles indi- 
cate the numerically calculated modes pi. In the same fig- 
ure, the time evolution of the distribution is shown for differ- 
ent times between t = 0.1ms and 0.6 s. For times greater 
than about 0.15 s, there is only very little change observed. 
As a function of time, the evolution of the single mode oc- 
cupations is shown in Fig. |4] We observe that the system 
very quickly, after about 5 ps, evolves to a quasistationary 
state, and that the subsequent drift to the equilibrium distri- 
bution takes roughly ten times longer In passing we note that 




f [s] 

FIG. 4: The normalized momentum-mode occupation numbers 
n{t;p)/niL, corresponding to those shown in Fig. [3] as functions 
of time. Shown are the populations of the modes with p = pi = 
2Ns/Lsm{iTr/Ns), i = 0, 1, iVs/2, and one has n{t;-p) = 
n{t;p). A fast short-time dephasing period is followed by a long 
quasistationary drift to the final equilibrium distribution. Notice the 
double-logarithmic scale. 



the mean-field Hartree-Fock (HF) approximation, for which 
= S'' = in Eqs. ( |45l l, conserves exactly all mode occu- 
pations and no equilibration is seen. 

In order to estimate to which extent the final distribution 
approaches that of the actual equilibrium state of the gas, we 
fitted the distribution to the Bose-Einstein-like form n{t; p) — 
[exp{(ci;(p) — /i)/fcB8(t;p)} — 1]^^, with a p-dependent tem- 
perature variable Q[t\p). Here was derived from the 
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FIG. 5: Momentum and time dependent temperature variable 
0{t;p) obtained by fitting the distribution n{t;p) — [exp{{uj{p) — 
li)/kBO{t;p)} — 1]^^ to the distribution obtained from the results 
shown in Fig. |4l for different, equally spaced times between t = 
and t = 0.6 s. One observes that, during the quasistationary pe- 
riod, 0.01 s< t < 0.1 s, no temperature can be associated to the 
distribution. Only at very large times, Q becomes approximately p- 
independent. 



time-derivatives of the statistical function F{t, t';p) att ~ t'. 
If a Bose-Einstein distribution is approached the temperature 
can be obtained from the slope of log(n^^ + 1) and the chem- 
ical potential /i from its value at uj = 0. Fig. |5] shows 
Q{t;p) for t = 0...0.6s. Obviously, during the quasistation- 
ary drift period, no temperature can be attributed to n{t;p), 
while, for large t, Q becomes approximately p-independent 
iFttIi . We deduce an approximate final temperature from 
e(0.6s; 128/L) = T = 0.35nK with ^ = LOSgiD^i for 
the above given values of (7id and rii, which, hence, deviates 
from the HFB result /i = giuni by 8% only. 

We furthermore studied, in the spirit of Refs. iItoI 
I72I1 . the time dependence of the ratio of the envelopes 
of the unequal-time correlation functions, specifically, 
at;p) = [(Fn(i,0;p)2 + Fi2(t,0;p)2)/(pii(t,0;p)2 + 
Pi2(i, 0;p)^)]^/^/n(i;p). Fig. |6]shows ^, for four different 
momentum modes, as a function of time. Due to the normal- 
ization with respect to n{t;p) all ^{t;p) are of the same order 
of magnitude. However, they show a distinct time evolution 
during the dephasing period, before they settle to a constant 
value during the quasistationary drift. ^ is a measure of the in- 
terdependence of the statistical and spectral functions, which, 
in thermal equilibrium, are connected through the fluctuation- 
dissipation relation [JOJ. In the momentum-frequency do- 
main, this relation is given in Eq. ( |79] l. Hence, the stationary 
^ indicates that F and p are linked to each other long before 
the momentum distribution becomes thermal. 

A further signature is found when comparing the kinetic 
and interaction contributions to the total energy as shown in 
Fig. |7] During the quasistationary drift, these contributions 
are constant and of the same order of magnitude, calling in 
mind the virial theorem. 

In summary, during the drift period, the system is not yet 
in equilibrium as far as the momentum distribution and tem- 
perature is concerned, but shows important characteristics of 
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FIG. 6: 



(Color online) Ratio of the envelopes of the 



unequal-time correlation functions, ^{t;p) = [(Fn (t, 0;p) + 
F^2it,0;pf)/{piiit,0;pf + p^2{t,Q-pf)f/^/n{t-p), for four 
different momentum modes, as a function of time. Due to the nor- 
malization with respect to n{t\ p) all ^(t; p) are of the same order of 
magnitude. ^ is a measure of the interdependence of the statistical 
and spectral functions, and its settling to a constant value during the 
quasistationary drift period indicates that these functions become, as 
in thermal equilibrium, connected through a fluctuation-dissipation 
relation. 
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FIG. 7: (Color online) Evolution of the kinetic and interaction con- 
tributions to the total energy of the gas. During the quasistationary 
drift these contributions assume the same order of magnitude, calling 
in mind the virial theorem. 



a system close to equilibrium. 



C. Evolution of the classical gas 

In SectionlUwe have discussed in detail the distinction be- 
tween the quantum and classical statistical contributions to the 
2PI effective action and to the dynamic equations. In the fol- 
lowing we compare the predictions for the classical statisti- 
cal theory with those presented in the preceding section, and 
pointing to the distinct differences. 

We solved Eqs. ( |45] l, with the self-energies now given by 
Eqs. ( |77] ). for the initial conditions given in Section lTlI Al first 
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FIG. 8: (Color online) The normalized momentum-mode occupation 
numbers n[t\p) /n\L for the classical gas (black solid lines) com- 
pared to their quantum counterparts from Fig. |4] (red dashed lines), 
as functions of time. Shown are the populations of the modes with 
p — Pi — 2Ns/ Lsiniin /Na), i = 0, 1, and one has 
n{t; — p) = n(t;p). In contrast to the quantum statistical evolution 
there is no quick dephasing in the classical case, such that the ini- 
tially empty modes become only subsequently filled, as is discussed 
in the main text. At large times, the classical and quantum gases 
necessarily evolve to different distributions. 



with a — 0. For comparison we also consider the case a = 1 
such that the classical and quantum initial correlation func- 
tions F and p are identical. In the classical case, the self- 
energies as well as the equations (173] ) determining the cou- 
pling functions Ip ^ are lacking certain terms compared to 
their quantum counterparts, cf. Eq. dTTj i. 

The time evolution of the initially Gaussian far-from- 
equilibrium momentum distribution of the classical gas is 
shown in Fig. [8] The mode occupations are shown as (black) 
solid lines, and for better comparison, we have quoted the 
quantum evolution from Fig. |4]as (red) dashed curves. One 
observes that the time evolution of the modes with occupa- 
tion number n{t;p) > 1, i.e., n{t;p)/niL > 10^^, is, for 
most of the time, identical to that obtained in the quantum 
case, confirming condition (fTSl ). As expected, for the strongly 
populated modes of a weakly interacting bosonic gas quan- 
tum fluctuations do not play a significant role for not too large 
times. Only when the evolution approaches the equilibrium 
state, the differences between quantum and classical statistics 
are expected to lead to a Bose-Einstein and classical distribu- 
tion, respectively. We point out that, although we have cho- 
sen, for our comparison, the same initial occupation numbers 
in the two cases, the total energies are different since the ini- 
tial correlation functions differ according to Eq. ( l83b . Hence, 
also the final-state occupation numbers of the low momentum 
modes can differ 

We finally point to the substantial differences in the short- 
time evolution of the weakly populated modes. While, in 
the quantum gas, the large-momentum mode populations are 
all growing at the same rate, the modes of the classical gas 
become populated much more gradually. The quantum-gas 
modes are occupied by 0.01 and 1 particle per mode already 
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FIG. 9: (Color online) The normalized momentum-mode occupation 
numbers n[t\p) /n\L for the classical gas (black solid lines) com- 
pared to their quantum counterparts from Fig. |4] (red dotted lines), 
as functions of time. All parameters are chosen as in Fig. [8] except 
for a — 6.5 • 10* m^^. The dashed (blue) lines show n(t,p) — 1/2, 
i.e., after subtracting the additional flat initial distribution which sim- 
ulates the quantum "zero point fluctuations". 



between 0.5 and 1 /is, while the classical modes need up to ten 
times longer. The distinct quantum behaviour of the modes 
can be understood as follows. At the energies present in 
such an ultracold gas, the atomic interactions are essentially 
pointlike, i.e., the range of the potential is not resolved and 
the coupling function or scattering amplitude is constant over 
the range of relevant momenta. Hence, in a single scattering 
event, the distance of two atoms is localized to zero, such that 
the relative momentum of the atoms is completely unknown 
immediately after the collision. This means that the transfer 
probability of the atoms is the same for any final momentum 
mode, and this is observed as a quick, collective population in 
the respective, so far essentially unoccupied modes. 

For comparison we repeated our calculations for a non- 
local interaction potential. In the momentum domain, this 
coiTesponds to a coupling function which is cut off at large 
momenta, and we chose the cutoff within the range of the mo- 
menta shown explicitly in Fig. [8] In this case we find that the 
quantum evolution is modified such that it becomes similar to 
that of the classical gas. In particular, all modes above the 
cutoff populate gradually one after each other 

The differences between the quantum and classical evo- 
lutions shown above depend, however, considerably on the 
choice of initial conditions. To compare the characteristics 
of the evolutions which are independent of the initial choice 
of F, we have repeated the classical calculations for an ini- 
tial momentum distribution where, as compared to before, a 
constant occupation number 1 /2 has been added. Hence, we 
chose a = 1 in the initial values of F, Eq. ( [83] l, as in the quan- 
tum case, such that F is identical for xq = i/q = to in the clas- 
sical and quantum cases. The results are shown in Fig.|9] The 
(red) dotted lines show, again, the quantum evolution, while 
the classical mode populations for the same initial conditions 
for F and p are shown as solid (black) lines. Subtracting 1/2 
from each Faa{t, t;p) gives the dashed (blue) lines. We find. 
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that during the initial period the evolution of the variation of 
the high-momentum modes with respect to their initial occu- 
pation is identical to the quantum evolution of the occupation 
numbers. At intermediate times, however, there are devia- 
tions which lead to occupation numbers up to a percent lower 
than 1/2. Although the chosen initial conditions which, in the 
classical case correspond to a base occupation of each mode 
with 1/2 atom, seem unphysical, our results show that for the 
dilute, weakly interacting gas under consideration, there are 
differences only in those modes which, in the mean, are popu- 
lated with less than one atom. Quantum statistical fluctuations 
play a role only for these modes. 

The situation is very different for an ultracold gas in the 
strongly interacting regime. We obtain this from the case dis- 
cussed so far by reducing the total number of atoms by a factor 
of 100 and increasing the dimensionless interaction parame- 
ter 7 by lO**. In this case the 'classicaUty' condition ( iTSl l is 
clearly violated and we expect strong corrections due to quan- 
tum fluctuations. We emphasize that here our nonperturba- 
tive approximation is crucial in order to be able to consider 
such strongly interacting systems. Note that the structure of 
the classical equations of motion, Eqs. ( l45l l. with the self- 
energies given by Eqs. dTTl i. is such that they are invariant un- 
der a simultaneous rescaling of 7 and ni which leaves 771^ un- 
changed. The quantum equations will, however change, with 
the specific quantum terms becoming more and more impor- 
tant with growing 7 and reduced density ni. We compare the 
evolution, again, for identical initial conditions as in the case 
shown in Fig. |9] The deviations between quantum and clas- 
sical statistical evolution are now quantitatively substantial. 
While many qualitative aspects are similar to those discussed 
above, now all momentum modes show clear deviations. As a 
characteristic example, we present, in Fig.[Tol the ratio of the 
envelopes of the unequal-time correlation functions, in anal- 
ogy to Fig. |6] for the quantum and classical evolution of the 
strongly interacting gas. One observes that the unequal-time 
correlation function decays more rapidly in the classical statis- 
tical case. Our results illustrate that the quantum system keeps 
much longer the information about the initial conditions, here 
roughly by an order of magnitude in evolution time. On the 
other hand, the example also shows that the classical statisti- 
cal system still behaves qualitatively similar to the quantum 
gas, despite the low densities and strong interactions. We 
would like to emphasize that clearly distinguishing quantum 
and classical statistical fluctuations will typically require high 
precision, both from the calculational as well as experimental 
point of view. 



IV. CONCLUSIONS 

In this article we have studied the non-equilibrium dynam- 
ics of an ultracold Bose gas and compared the theoretical pre- 
dictions from a full quantum approach with those where only 
classical statistical fluctuations were taken into account. On 
the basis of functional techniques in quantum field theory we 
reviewed the difference in the generating functionals, the 2PI 
effective actions which determine the dynamic equations, be- 




FIG. 10: (Color online) Ratio of the envelopes of the 
unequal-time correlation functions, ^{t;p) = [(Fn (t, 0; p)^ + 
F^2it,0;pf)/{piiit,0;pf + pi2it,0;pf)Y^yn{t;p), for four 
different momentum modes, as a function of time. This figure is 
analogous to Fig. (6) but shows the difference of the quantum (black 
lines) and classical (blue lines) evolutions for a gas in the strongly 
correlated regime 7 = 15, for ni = 10^ m~^. The results show that 
in this regime, where quantum and classical statistical evolution dif- 
fer considerably, the quantum evolution conserves information about 
the initial conditions much longer, here roughly by an order of mag- 
nitude in time. 



tween the quantum and the classical cases. The functional 
descriptions for both cases show to be very similar In partic- 
ular, the crucial difference is the absence, in the classical ver- 
sus the quantum case, of certain coupling terms or vertices in 
the action which defines the theory and which enters the gen- 
erating functional of correlation functions. As a consequence, 
the classical generating functional is characterized by an im- 
portant reparametrization property, such that for the classical 
dynamics the effects of a larger self-interaction can always 
be compensated by a smaller density. Quantum corrections 
violate this invariance property. They become of increasing 
importance with growing scattering length or reduced density. 
We have used this to derive a 'classicality' condition, which is 
not based on thermal equilibrium assumptions and holds also 
for far-from-equilibrium dynamics. In particular, it is inde- 
pendent of the value of the macroscopic field (p. 

To illustrate the possibilities one has at hand, we studied 
the equilibration dynamics of a homogeneous ultracold gas of 
interacting sodium atoms in one spatial dimension. The gas 
is assumed to be, initially, in a state far from thermal equi- 
librium, characterized by a Gaussian momentum distribution 
centred at zero momentum. Extending on earlier results [3^ 
we find that the evolution takes place via two distinct periods, 
with a fast initial dephasing followed by a slow quasistation- 
ary drift to the final equilibrium distribution. While the gas, 
during the quasistationary drift, already shows important char- 
acteristics of a near-equilibrium situation like a fluctuation- 
dissipation relation and constancy of the total kinetic energy, 
no temperature can yet be associated to the momentum distri- 
bution. 

We compared this evolution to the evolution of a classi- 
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cal gas with different initial conditions. If the initial mode 
occupation numbers are chosen as in the quantum case, no 
short-time dephasing period is found. The modes are rather 
populated subsequently, with the quasistationary drift setting 
in gradually. The long-time approach to the final equilibrium 
state expectedly shows a distinct behaviour compared to the 
quantum gas. We moreover showed that the dephasing in 
the short-time evolution of the quantum gas can be approxi- 
mately simulated by a classical gas if one chooses appropri- 
ate initial conditions. Our results demonstrate explicitly that 
quantum fluctuations only play a role for modes whose oc- 
cupation number is sufficiently small. We found that the dy- 
namics of the weakly interacting one-dimensional gas is al- 
most purely classical. In contrast, for a strongly interacting 
one-dimensional gas we observed substantial quantitative de- 
viations. There are distinctive properties of quantum versus 
classical statistical dynamics, as e.g. that, given identical ini- 
tial values for the correlation functions, in the quantum evo- 
lution, information about the initial conditions is conserved 
much longer as it would be expected from a calculation in 
classical approximation. Such differences may be expected 
to be of interest in the context of strongly correlated ultra- 
cold atomic gases, e.g. in lower-dimensional traps or near 
Feshbach scattering resonances. However, we emphasize that 
high precision measurements together with accurate calcula- 
tions are typically required to be able to clearly distinguish be- 
tween effects of quantum and classical statistical fluctuations. 
We think that the nonperturbative methods presented here can 
be a very valuable tool to precisely identify the effect of quan- 
tum fluctuations on the time evolution of ultracold quantum 
gases observed in present-day and near-future experiments. 
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APPENDIX A: INITIAL CONDITIONS 

In this appendix we provide explicit expressions for a gen- 
eral Gaussian initial-state density matrix po [(^o + <po/2, fvi — 
(po/2] and a Gaussian probability functional W^[(y9o,7ro] en- 
tering the generating functionals (fTSl l and (|58] l, respectively. 
These specify the initial conditions for the dynamic equations 
(l45T l. (ISOl l. We provide expressions for a spatially homoge- 
neous system, as considered in Sect. |III] 

The most general Gaussian initial density matrix takes, in 



the representation (fT2] l of the fields, the form 

PD[<y5o + '^o/2,Vo - <^o/2] = r, ]- i- exp I 

.■qa 



y^o,a<yab'Po,b 



1 Va 



(Al) 



where summation over a and b is implied in the exponent. 
Since pu involves a factor of the above form above for each 
momentum mode, mode indices at the fields and parameters, 
as well as summation over momenta in the exponent have been 
neglected. In order to reflect the symmetry ^ between the 
field and the canonical momentum, the six parameters £^a, Va, 
and (Ta are reduced to three independent parameters through 
the conditions 



ai = a2 — (J, 

t2 _ „2 I 2 



m = -12^2/ £.1- 



(A2) 



Using the definition of the initial statistical correlation func- 
tions in terms of the initial density matrix. 



Fab{0,0) = pi3(to){$a(0),$f,(0)} 



PO,a(PO,b, 

(A3) 



and inserting unit operators J[d93]|(^)(</?| and/or J [dvr] |7r) (7r| 
one finds that the three free parameters are related to the initial 
correlation functions as follows 



£,1 — Fiiih, h) 
£im = Fi2{to,tQ) 



t>0,l4'0,2, 
t>0,2'p0,2- 



(A4) 



Again, all mode labels have been suppressed. The initial spec- 
tral functions Pab{0, 0) do not enter the density matrix as they 
are fixed, in the quantum and classical statistical cases, by the 
commutation relations and the Poisson brackets, respectively. 

Note that the initial conditions chosen in Sect. lIII Al for the 
numerical evaluation of the dynamics of a one-dimensional 
Bose gas, with Fii(0,0) = i^22(0,0), i^i2(0,0) = and 
(f,n^a = con-espond to £f ^ ^| = - Fn{0, 0), rja = 0, 
CT = 2^2, such that 



Poifa + <^o/2, ifa - <Po/2] 
1 



■ exp 



2 

27re2 i 2^2 ^0,a 2 



1 p 

^2 S ~2 



(A5) 



We close this appendix by providing the expression for the 
probability functional VF[v3oj ""o] obtained through the inverse 
of the Fourier transform ( |59] l: 



W[ipo,itq] 
1 



-^0,a — TTo.a 



2er 



(A6) 
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(Summation over a in the exponent). It is straightforward to with \Fab\ = Fii{0, 0)^22(0, 0) — -Fi2(0, 0)^21 (0, 0). 
show that the above expression is symmetric under exchange 
of (^o,a and 7ro,a, and that, by inserting 7ro,a = ia'^^ipo^, one 
obtains the expression 

+ i^22(0, 0) - 2i^i2(0, 0) ¥'0,1^0,2] }, (A7) 
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